The Argyris isogeometric space on unstructured multi-patch planar domains

Multi-patch spline parametrizations are used in geometric design and isogeometric analysis to represent complex domains. We deal with a particular class of $C^0$ planar multi-patch spline parametrizations called analysis-suitable $G^1$ (AS-$G^{1}$) multi-patch parametrizations (Collin, Sangalli, Takacs; CAGD, 2016). This class of parametrizations has to satisfy specific geometric continuity constraints, and is of importance since it allows to construct, on the multi-patch domain, $C^1$ isogeometric spaces with optimal approximation properties. It was demonstrated in (Kapl, Sangalli, Takacs; CAD, 2018) that AS-$G^1$ multi-patch parametrizations are suitable for modeling complex planar multi-patch domains. In this work, we construct a basis, and an associated dual basis, for a specific $C^1$ isogeometric spline space $\mathcal{W}$ over a given AS-$G^1$ multi-patch parametrization. We call the space $\mathcal{W}$ the Argyris isogeometric space, since it is $C^1$ across interfaces and $C^2$ at all vertices and generalizes the idea of Argyris finite elements to tensor-product splines. The considered space $\mathcal{W}$ is a subspace of the entire $C^1$ isogeometric space $\mathcal{V}^{1}$, which maintains the reproduction properties of traces and normal derivatives along the interfaces. Moreover, it reproduces all derivatives up to second order at the vertices. In contrast to $\mathcal{V}^{1}$, the dimension of $\mathcal{W}$ does not depend on the domain parametrization, and $\mathcal{W}$ admits a basis and dual basis which possess a simple explicit representation and local support. We conclude the paper with some numerical experiments, which exhibit the optimal approximation order of the Argyris isogeometric space $\mathcal{W}$ and demonstrate the applicability of our approach for isogeometric analysis

Construction of analysis-suitable G1G^1 planar multi-patch parameterizations

Isogeometric analysis allows to define shape functions of global $C^{1}$ continuity (or of higher continuity) over multi-patch geometries. The construction of such $C^{1}$-smooth isogeometric functions is a non-trivial task and requires particular multi-patch parameterizations, so-called analysis-suitable $G^{1}$ (in short, AS-$G^{1}$) parameterizations, to ensure that the resulting $C^{1}$ isogeometric spaces possess optimal approximation properties, cf. [7]. In this work, we show through examples that it is possible to construct AS-$G^{1}$ multi-patch parameterizations of planar domains, given their boundary. More precisely, given a generic multi-patch geometry, we generate an AS-$G^{1}$ multi-patch parameterization possessing the same boundary, the same vertices and the same first derivatives at the vertices, and which is as close as possible to this initial geometry. Our algorithm is based on a quadratic optimization problem with linear side constraints. Numerical tests also confirm that $C^{1}$ isogeometric spaces over AS-$G^{1}$ multi-patch parameterized domains converge optimally under mesh refinement, while for generic parameterizations the convergence order is severely reduced

A family of C1C^1 quadrilateral finite elements

We present a novel family of $C^1$ quadrilateral finite elements, which define global $C^1$ spaces over a general quadrilateral mesh with vertices of arbitrary valency. The elements extend the construction by (Brenner and Sung, J. Sci. Comput., 2005), which is based on polynomial elements of tensor-product degree $p\geq 6$, to all degrees $p \geq 3$. Thus, we call the family of $C^1$ finite elements Brenner-Sung quadrilaterals. The proposed $C^1$ quadrilateral can be seen as a special case of the Argyris isogeometric element of (Kapl, Sangalli and Takacs, CAGD, 2019). The quadrilateral elements possess similar degrees of freedom as the classical Argyris triangles. Just as for the Argyris triangle, we additionally impose $C^2$ continuity at the vertices. In this paper we focus on the lower degree cases, that may be desirable for their lower computational cost and better conditioning of the basis: We consider indeed the polynomial quadrilateral of (bi-)degree~$5$, and the polynomial degrees $p=3$ and $p=4$ by employing a splitting into $3\times3$ or $2\times2$ polynomial pieces, respectively. The proposed elements reproduce polynomials of total degree $p$. We show that the space provides optimal approximation order. Due to the interpolation properties, the error bounds are local on each element. In addition, we describe the construction of a simple, local basis and give for $p\in\{3,4,5\}$ explicit formulas for the B\'{e}zier or B-spline coefficients of the basis functions. Numerical experiments by solving the biharmonic equation demonstrate the potential of the proposed $C^1$ quadrilateral finite element for the numerical analysis of fourth order problems, also indicating that (for $p=5$) the proposed element performs comparable or in general even better than the Argyris triangle with respect to the number of degrees of freedom

Adaptive isogeometric methods with C1C^1 (truncated) hierarchical splines on planar multi-patch domains

Isogeometric analysis is a powerful paradigm which exploits the high smoothness of splines for the numerical solution of high order partial differential equations. However, the tensor-product structure of standard multivariate B-spline models is not well suited for the representation of complex geometries, and to maintain high continuity on general domains special constructions on multi-patch geometries must be used. In this paper we focus on adaptive isogeometric methods with hierarchical splines, and extend the construction of $C^1$ isogeometric spline spaces on multi-patch planar domains to the hierarchical setting. We introduce a new abstract framework for the definition of hierarchical splines, which replaces the hypothesis of local linear independence for the basis of each level by a weaker assumption. We also develop a refinement algorithm that guarantees that the assumption is fulfilled by $C^1$ splines on certain suitably graded hierarchical multi-patch mesh configurations, and prove that it has linear complexity. The performance of the adaptive method is tested by solving the Poisson and the biharmonic problems

A locally based construction of analysis-suitable G1G^1 multi-patch spline surfaces

Analysis-suitable $G^1$ (AS-$G^1$) multi-patch spline surfaces [4] are particular $G^1$-smooth multi-patch spline surfaces, which are needed to ensure the construction of $C^1$-smooth multi-patch spline spaces with optimal polynomial reproduction properties [16]. We present a novel local approach for the design of AS-$G^1$ multi-patch spline surfaces, which is based on the use of Lagrange multipliers. The presented method is simple and generates an AS-$G^1$ multi-patch spline surface by approximating a given $G^1$-smooth but non-AS-$G^1$ multi-patch surface. Several numerical examples demonstrate the potential of the proposed technique for the construction of AS-$G^1$ multi-patch spline surfaces and show that these surfaces are especially suited for applications in isogeometric analysis by solving the biharmonic problem, a particular fourth order partial differential equation, over them

C1C^1-smooth isogeometric spline functions of general degree over planar mixed meshes: The case of two quadratic mesh elements

Splines over triangulations and splines over quadrangulations (tensor product splines) are two common ways to extend bivariate polynomials to splines. However, combination of both approaches leads to splines defined over mixed triangle and quadrilateral meshes using the isogeometric approach. Mixed meshes are especially useful for representing complicated geometries obtained e.g. from trimming. As (bi)-linearly parameterized mesh elements are not flexible enough to cover smooth domains, we focus in this work on the case of planar mixed meshes parameterized by (bi)-quadratic geometry mappings. In particular we study in detail the space of $C^1$-smooth isogeometric spline functions of general polynomial degree over two such mixed mesh elements. We present the theoretical framework to analyze the smoothness conditions over the common interface for all possible configurations of mesh elements. This comprises the investigation of the dimension as well as the construction of a basis of the corresponding $C^1$-smooth isogeometric spline space over the domain described by two elements. Several examples of interest are presented in detail